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Modeling  and  Control  of  Combustion  Instability  Using  Fuel 

Injection  1 

J.P.  Hathout,  A.M.  Annaswamy  and  A.F.  Ghoniem 
Department  of  Mechanical  Engineering 
Massachusetts  Institute  of  Technology 
Cambridge,  MA  02139,  USA 


Abstract 


Active  control  using  periodic  fuel  injection  has  the  po- 
tenlial  of  suppressing  combusLion  instability  without 
radically  changing  the  engine  design  or  sacrificing  per- 
formance. In  this  paper,  we  carry  out  a study  of  opti- 
mal model-based  control  of  combustion  instability  us- 
ing fuel  injection.  The  model  developed  is  physically 
based  and  includes  the  acoustics,  the  heat-release  dy- 
namics, their  coupling,  and  the  injection  dynamics.  A 
heat-release  model  with  fluctuations  in  the  flame  sur- 
face area  as  well  as  in  the  equivalence  ratio  is  derived. 
We  show  that  area  fluctuations  coupled  with  the  veloc- 
ity fluctuations  drive  longitudinal  modes  to  resonance 
caused  by  phase-lag  dynamics,  while  equivalence  ratio 
fluctuations  can  destabilize  both  longitudinal  and  bulk 
modes  caused  by  timc-dclav  dynamics,  similar  to  ex- 
perimental observations.  The  dynamics  of  proportional 
and  two-position  (on-off)  fuel  injectors  are  included  in 
the  model.  Using  the  overall  model,  tw'O  different  con- 
trol designs  are  proposed.  The  first  is  an  LQG/LTR 
controller  w'here  the  time-delay  effect  is  ignored,  and 
the  second  is  a Posi-Cast  controller  w'hich  explicitly  ac- 
counts for  the  delay.  Injection  at  (i)  the  burning  zone 
and  (ii)  further  upstream  is  considered.  The  character- 
istics of  fuel  injectors  including  bandwidth,  authority 
(pulscd-fuel  flow  rate),  and  whether  it  applies  a propor- 
tional ora  two-position  (on-off) injection  are  discussed. 
We  show  that  increasing  authority  and  bandwidth  result 
in  improved  performance.  Injection  at  (ii)  compared  to 
(i)  results  in  a trade-off  between  improved  mixing  and 
increased  time-delay.  We  also  note  that  proportional 
injection  is  more  successful  than  on-off  injection  since 
the  former  can  modulate  both  amplitude  and  phase  of 
the  control  fuel. 


1This  work  is  sponsored  in  part  by  the  National  Science  Founda- 
tion. contract  no.  PCS  9713415.  and  in  part  by  the  Office  of  Naval 
Research,  contract  no.  N00014-99-1-0448. 


1 Introduction 

Active  control  using  periodic  fuel  injection  has  been 
recognized  as  a promising  technology  to  abaLe  combus- 
tion instability  in  practical  systems  [l]-[6].  Most  cur- 
rent designs  are  based  on  phase-shift  algorithms  which, 
despite  their  success,  arc  not  optimal  in  terms  of  fuel 
consumption,  settling  time,  and  robustness.  Phase-shift 
controllers,  due  to  their  limited  dynamics,  are  known  to 
succeed  only  within  a limited  frequency  band  (around 
the  unstable  frequency),  and  over  a small  range  of  op- 
erating conditions.  In  this  paper,  we  present  alterna- 
tive control  designs  using  LQG/LTR  and  Posi-Cast  con- 
trol methods  which  are  based  on  a physical  model,  and 
demonstrate  their  advantages. 

Control  algorithms  must  be  based  on  an  accurate  de- 
scription of  the  system  dynamics  to  achieve  optimal 
performance.  Much  effort  has  gone  into  modeling 
acoustics  [ l],[7]-[  1 0],  Modeling  of  heat-release  dy- 
namics and  its  coupling  with  acoustics,  mixing  dynam- 
ics, and  actuator  dynamics,  despite  their  crucial  role  in 
the  instability  mechanism,  have  not  been  dealt  with  to 
the  same  extent.  Recognizing  that  heat-release  oscilla- 
tions can  be  caused  by  fluctuations  of  the  flame  surface 
area  due  to  velocity  perturbation  or  the  equivalence  ra- 
tio due  to  pressure  perturbation,  we  derive  a model  that 
incorporates  both  for  the  case  of  a weakly  turbulent  pre- 
mixed flame.  We  show  how  certain  acoustic  modes  can 
couple  with  either  or  both  forms  of  heat  oscillations. 

The  feasibility  of  control  design  is  tightly  correlated 
w'ith  the  performance  of  sensors  and  actuators.  While 
high-bandwidth  devices  of  the  first  are  available,  c.g., 
pressure  transducers  and  heat-release  sensors,  actuation 
by  means  of  fuel  injection  is  hindered  by  low  band- 
width, limited  authority,  and  nonlinearities  (in  the  form 
of  dead-zone,  saturation  and  on-off  effects)  [11],  The 
impact  of  ideal  actuators  was  studied  extensively  in 
[12],  The  ability  of  the  actuator  to  convey  the  control 
signal  faithfully  to  the  combustor,  thus  introducing  the 
correct  forcing,  is  key  to  stable  performance.  This  how- 
ever is  often  not  the  case  due  to  (i)  non-ideal  injection 


Paper  presented  at  the  RTO  AVT  Symposium  on  “Active  Control  Technology  for 
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dynamics,  and  (ii)  injection  locations  which  may  im- 
pose a delay  in  the  control  input.  To  capture  the  contri- 
bution of  these  effects  in  the  control  design,  we  develop 
a model  of  the  injector  dynamics  and  incorporate  it  in 
the  overall  analysis  of  the  system. 

In  Sec.  2,  we  model  the  heat  release  dynamics  forced  by 
velocity  and  equivalence  ratio  fluctuations,  the  acoustic 
dynamics,  the  coupling  and  mixing  dynamics,  their  ef- 
fect on  stability,  and  a generic  model  of  typical  solenoid 
injectors  which  captures  its  relevant  features.  In  Sec. 
3,  control  is  developed  for  a multi-mode  combustor 
and  injection  at  as  well  as  upstream  of  the  burning 
zone  is  implemented.  For  the  former,  an  optimal 
Linear-Quadratic  Gaussian  with  Loop-Transfer  Recov- 
ery (LQG/LTR)  controller  is  used,  and  for  the  latter,  a 
novel  technique  known  as  Posi-Cast  control  that  accom- 
modates combustion  systems  with  delays  is  developed 
and  implemented. 


2 Combustor  Model 

Physically-based  modeling  has  been  shown  to  be  indis- 
pensable for  optimal  control  design  for  suppression  of 
combustion  instability  [13].  The  dynamics  leading  to 
combustion  instability  is  known  to  occur  due  to  the  res- 
onant coupling  between  acoustics  and  heat  release.  This 
coupling  may  occur  through  flame-area  fluctuations,  as 
shown  experimentally  in  [1],  and  by  modeling  in  [14]. 
It  can  also  occur  through  mixture  inhomogeneity  (reac- 
tants’ equivalence-ratio  fluctuations),  observed  experi- 
mentally in  [1 5]-[l 7],  and  modeled  in  [18,  19],  While 
both  mechanisms  can  destabilize  longitudinal  modes 
[4,  7,  16,  20].  bulk  modes  are  more  strongly  affected 
by  the  latter  [3.  18].  A rigorous  model  for  the  response 
of  heat-release  dynamics  to  these  oscillations  is  shown 
next.  This  is  followed  by  a description  of  acoustics, 
coupling  mechanisms,  and  fuel-injector  dynamics. 

2.1  Heat-Release  Dynamics 

The  response  of  the  heat -release  rate  to  perturbations  in 
the  flow  velocity,  u,  as  well  the  equivalence  ratio,  (f>  are 
modeled  in  this  section.  The  flame  kinematics  equation 
is  used  to  derive  this  model.  The  following  assump- 
tions are  made:  (i)  The  flame  is  a thin  interface  sepa- 
rating reactants  and  products  and  is  insensitive  to  pres- 
sure perturbations  [1],  The  flame  can  model  turbulent 
premixed  combustion  if  conditions  of  high  Damkohler 
number  and  weak  to  moderate  turbulence  intensity  pre- 
vail [21]-[23].  (ii)  The  flame  is  w'eakly  convoluted,  i.e., 
it  can  be  described  by  a single-valued  function,  £(r,  t), 
describing  the  instantaneous  location  of  its  surface. 

Under  the  above  assumptions,  the  heat-release  model  is 


of  the  form  1 


where  Q is  the  heat-release.  Su  is  the  burning  velocity, 
k{4>)  = 2icpuSu{<j>)\hr(<l>).  pu  is  the  density  of  the 
unburnt  mixture,  and  A/i  ,.  is  the  heat  of  reaction. 


Assuming  negligible  velocity  component  in  the  radial 
direction,  and  linearizing  around  nominal  values  u,  S u. 
and  £(>■),  denoting  (.)  and  (.)'  as  steady  and  perturba- 
tion, respectively,  we  get 


0? 

Ot 


tj  W + OidSu 
" Or  Or  dtp 


wdth  boundary  conditions" 


(3) 


?{K,t)  = 0 Vf, 


£V,U)  = U Vc, 


while 


Q'(t) 


k / f'( r,t)dr  | di4> , (4) 

J o 


where3 


di  = 2ttpu 


k = 2w/juSuAhr, 


ll  is  worth  noting  that  the  flame  area  fluctuation.  .4 U is 

given  by  A'f(t)  — 2 tt  fj*  £'(r,t)dr.  This  with  Eq.  (3) 
shows  that  the  flame  area  is  affected  by  both  ur  and  ©', 
and  the  area  in  turn  impacts  Qr  as  shown  in  Eq.  (4). 
This  also  shows  that  tf>'  affects  Q'  directly,  and  indi- 
rectly through  the  area  fluctuations. 


Equation  (3)  can  be  manipulated  further  and  solved  for 
in  the  Laplace  domain  as: 


C'(r,«)= 


u'{$) 


0£  dSu 
Or  dtp 


-i  R-r) 


(5) 


where  s is  the  Laplace  operator.  Differentiating  Eq.  (4) 
with  respect  to  time,  using  Eq.  (3),  and  integrating  over 
[0,  R],  we  obtain  that 


Q =7c  .Ru' - S„£'( 0,t) 


r 21  21a 

Jo  dr  dd> 


(pr dn-\di(j> . (6) 


1 We  consider  here  a flame  stabilized  over  a perforated  plate,  Ilh 
the  radius  of  the  perforation. 

2lt  should  be  noted  that  with  the  appropriate  change  in  coordinates 
and  boundary  conditions,  Eq.  (3>  can  also  represent  flames  stabilized 
behind  a gutter  [24],  or  a dump  [22], 

3The  factor  1_  is  positive  and  |—  is  also  positive  when 
0 < 1. 
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Taking  the  inverse  Laplace  of  Eq.  (5)  at  r = 0,  and 
substituting  in  Eq.(6),  after  some  manipulations,  we  get 


and  Tf  — R/Su  is  the  characteristic  propagation  delay 
of  the  flame  surface  into  the  reactants  flow.  Note  that  for 
the  class  of  flames  considered  in  the  paper,  the  slope  at 
the  flame  tip,  which  is  typically  conical,  is  zero,  there- 
fore the  third  term  in  the  RHS  of  Eq.  (7)  can  be  omitted. 

The  effect  of  u'  on  Q'  was  previously  considered  in 
[14]  where  the  time-delay  effects  were  assumed  to  be 
small.  The  effect  of  q>'  on  ()'  was  addressed  in  [18] 
where  a lumped  parameter  approach  was  used  where 
once  again,  the  time-delay  effects  were  neglected.  In 
contrast,  we  have  analyzed  the  effect  of  both  pertur- 
bations on  the  flame  dynamics  u'here  the  impact  of  a 
time-delay,  Tf,  is  taken  into  account. 


2.2  Acoustics 

Acoustic  resonance  of  a combustor  can  be  captured  by 
analyzing  the  wave  equation  and  linearizing  the  conser- 
vation equations  [20] 


d~p'  d~p' 

W~crd^ 


(y  -l)q'(z,t),  (8) 


where  p is  the  pressure,  and  c is  the  mean  speed  of 
sound.  We  consider  flames  localized  close  to  the  an- 
choring plane,  q'(x,t)  — q'{t)5{x  — Xf). 


Using  an  expansion  in  basis  functions 


p'(x,t)  = p ^2  ti  (x)  7li  (0  * (9) 

i— 1 


where  in  most  cases,  ify.  (x)  = sin  (k;x  + <pi0),  hi  and 
d>io  determined  from  the  boundary  conditions,  and  per- 
forming a weighted  spatial  averaging,  the  modal  ampli- 
tudes can  be  shown  Lo  follow  [20]: 


2.3  Coupling  Dynamics 

Perturbations  in  the  flame  area  and  the  equivalence  ratio 
w’hich  couple  in  a resonant  way  with  the  acoustic  field 
are  caused  by  velocity  or  pressure  oscillations.  Pertur- 
bations in  the  velocity  and  pressure  are  related  to  the 
perturbation  in  acoustic  field  via  the  energy  equation 
[20] 


dp' 


_du' 


Ht  + ^~dx 


= (7  -IV. 


(11) 


which  can  be  integrated  over  the  combustor  length  to 
obtain  (using  Eq.  (9)) 


n 

u'  = ^C,:  T)i  +ea„q\  (12) 

i=  1 


where  -£<  (a '/),  u0  — (y  - l)/y p,  and  0 rep- 

resents the  effect  of  the  velocity  ahead  and  behind  the 
flame  on  q' . This  equation  can  be  used  w'ith  Eq.  (7)  to 
couple  the  heat  release  to  the  velocity  perturbations. 

Perturbations  in  the  equivalence  ratio  at  the  fuel  source 
due  to  coupling  with  the  acoustics  can  be  modeled  start- 
ing from  the  mass  conservation  at  the  mixing  section 
which  results  in 


<f> 


(13) 


for  fluctuations  in  the  air  flow  only  (fuel  inlet  is  choked) 
or  in  the  fuel  flow  only  (air  inlet  is  choked).  Now  using 
the  momentum  equation: 


0, 


(14) 


and  //  as  in  Eq.  (9),  after  some  manipulations,  we  get 


p<i 


dii';(z) 

~ <L~'lh 


(15) 


according  to  fluctuations  in  the  air  flow'  only  or  in  the 
fuel  flow  only.  The  equivalence  ratio  at  the  burning 
zone  is  related  to  that  at  the  fuel  source  by  a convec- 
tive delay,  rs  = Ls/u,  as 

4>' = 4>'s(t  - tv),  (16) 


n 

Vi  - 2 Tjj  +U>i7fi  = ^Tbi  (/f,  ( 1 0) 

i= 1 

where  = ya0ibi(xf)/E,  E — tj’f  (x)  ilx,  7 is 
the  specific  ratio,  a„  — Tyf,  f represents  the  passive 
damping  ratio  in  the  combustor4,  L is  its  length,  and 

tOj  — kifi. 

4 Dissipation  in  a combustor  can  be  caused  by  heat  losses  in  the 
flame  zone  and  friction  due  to  viscous  effects. 


w'here  rn  = ts  — TCOmb , and  the  latter  is  the  combustion 
time  delay.  In  most  eases  rs  >>  reomi„  and  hcncc 
rn  ~ ts  is  an  acceptable  approximation.  We  should 
also  mention  that  in  this  analysis  the  radial  distribution 
of  the  local  equivalence  ratio  is  assumed  to  be  constant. 

It  is  worth  noting  also  that  instability  can  arise  in  com- 
bustors where  flue-gas  recirculation  is  used  lo  reduce 
NOx  formation,  e.g.  [25].  The  recirculated  products 
may  couple  with  the  acoustics  and  cause  perturbations 
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in  the  equivalence  ratio  which  can  be  described  in  a 
similar  manner  as  in  Eqs.  (15)and  (16).  Also,  when  in- 
jection is  used  as  actuation,  it  introduces  an  additional 
term  to  appear  in  the  RHS  of  Eq.  (16)  and  is  discussed 
in  Sec.  3.1. 

2.4  The  Overall  Model 

Combining  the  acoustics,  heat-release  dynamics,  and 
convective-lag  effects  in  Eqs.  (10),  (7)  and  (16),  re- 
spectively, we  obtain  the  following  equations: 

i}i  4-  2 Qu>t  r/i  +ujf  rij=-J-K{Ru'  — Su  u'(t)cIt 

A,.  Jt-Tf 

+d}<P's(t  - rs)}  -k/i  4>[.  ( t - rs),  (17) 

where  Ac  is  the  cross-sectional  area  of  the  combustor. 
Equation  (17)  shows  that  two  different  time  delays,  r r 
and  rs,  can  lead  to  instability,  one  arising  from  flame 
propagation  effects,  and  the  other  from  fuel  convection. 
The  RHS  of  Eq.  (17)  can  be  simplified  further  noting 
that  for  a certain  class  of  flames  r/  is  small  when  Su 
is  large  (e.g.,  in  a turbulent  flow),  and  when  the  char- 
acteristic size  of  the  flame,  R,  is  small.  Also,  around 
the  unstable  acoustic  frequency,  <p  scales  as  u>|0|,  and 
0(7Zds)  ~ O(di),  thus  the  second  and  third  terms  in 
the  RHS  can  be  neglected: 


in  [4,  26]  where  the  reactants  inlet  were  choked  im- 
plying no  <p ' fluctuations.  This  condition  arises  wrhen 
the  flame  is  located  in  a position  where  the  modal  am- 
plitude of  the  pressure  and  the  velocity  (represented 
by  b and  c)  are  of  the  same  sign.  Eor  example,  in  a 
closed-open  combustor  with  fixed  equivalence  ratio,  a 
quarter-wave  mode  is  always  stable  whereas  a three- 
quarter-wave  mode  gets  destabilized  when  the  flame  is 
located  between  the  node  and  the  anti-node  of  the  pres- 
sure mode  [20], 

If  on  the  other  hand,  > 0,  even  though  this  implies 
that  the  damping  is  positive,  perturbations  in  the  equiv- 
alence ratio  can  destabilize  the  combustor  due  to  the 
convective  time-delay  t„.  When  r,  = 0,  the  system  is 
stable.  As  rs  increases,  the  system  becomes  unstable  if 
r„  lies  between  t~  and  t+,  where 


\Jl-2Q±  \/T|/w4  - 4£?(1  - gJ 

t IV^  ) 

if 

2Cc\/i-C-  < < 1. 


i);  + 2Cw;  r/f  =-^-  + d\4>s[t  - r,)j  . (18) 

The  overall  model  of  the  combustor  is  represented  by 
Eqs.  (18),  (12).  and  (15)  when  multiple  modes  are 
present.  For  ease  of  exposition,  we  restrict  our  discus- 
sion to  the  case  when  only  one  mode  is  present,  a sim- 
ilar approach  can  be  extended  to  the  multi-mode  case. 
Substituting  for  the  coupling  dynamics  in  Eqs.  (12), 
(15)  and  (16),  and  assuming  9 = 0,  we  get 

rj  + 2Qu)ij  + u>~ri  = IT// + rV//(f  — r6.),  (19) 
where  T1  = %kR  and  T,  - ± 

Defining  the  equivalent  damping  Co  as 

_ 2Cu>  — r i 

^ ~ 2u! 

when  perturbations  in  the  flame  area  change  the  equiv- 
alent damping  in  the  combustor,  such  that  Co  < 0,  the 
oscillator  is  always  unstable.  In  this  case,  T i repre- 
sents a “negative”  damping  added  by  the  flame  area 
fluctuations^,  and  is  suspected  to  cause  the  instability 

-’The  notion  of  “negative"  and  “positive”  damping  is  used  here  to 
indicate  heat-release  perturbation  which  is  out-of-phase  and  in-phase 
with  the  pressure  oscillation,  respectively.  We  have  shown  in  previous 
work  [12]  that  the  effect  of  flame- area  perturbation  can  be  modeled  as 
a negative  damping  term.  Heat  release  can  be  modulated  using  active 
forcing  to  produce  a positive  damping  term,  as  shown  in  [19]. 


if|r2|  > u.'2,  there  exists  only  one  switch  to  instability' 
at  t~,  though  in  most  practical  cases,  we  have  ] T ^ | < 
w2  [18]. 

W'hen  Co  ~ 0,  the  above  instability  conditions  simplify 
to  [18,  19]:  Defining  the  acoustic  time  constant  ra  — 
2k/<jJ.  the  system  is  unstable  if 


n \ 1 rs  ^ n I 2 

' Ta  ' 


(20) 


n — 0,2,4....  if  To  < ur.  and  unstable  if  £*■  > b,  if 

r2  > w2. 

The  above  discussion  indicates  that  the  instability  pat- 
tern depends  on  the  distance  between  the  source  of  </>'., 
the  convective  velocity,  and  the  acoustic  frequency  if 
perturbations  in  the  equivalence  ratio  instigate  the  in- 
stability which  appears  to  the  instability  mechanism  in 
the  rigs  in  [3,  16,  17].  The  dependence  of  the  stability 
bands  on  r5/ra  can  be  judiciously  exploited  for  opti- 
mizing the  active-control  designs  [19]. 


2.5  Fuel  Injector  Dynamics 

2.5.1  A Proportional  Injector:  The  pulsating 

fuel  injector  delivers  oscillations  in  the  mass-flow  rate 
in  response  to  a voltage  input.  The  injector  system  con- 
sists of  an  electro-mechanical  part  and  a fluidic  part, 
where  in  the  former,  the  input  voltage  generates  an 
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Figure  1:  A schematic  of  a typical  injector. 


electro-magnetic  field  that  causes  a poppet  to  move 
against  a spring  (as  seen  in  Fig.  1).  The  motion  of 
the  poppet  controls  the  aperture  of  the  injector  allowing 
fluid  to  flow. 


The  electro-mechanical  part  relates  the  voltage  E to 
the  poppet  position  x through  the  electrical,  electro- 
magnetic, and  mechanical  components,  which  can  be 
modeled  as 


E 


m 


d~x 

dtt- 


v 

Fm 

*t*  kx 


iRe  + T V,  (21) 

Bel~,  (22) 

at 

Beli,  (23) 

Fm , (24) 


where  i is  the  current,  Fm  is  the  magnetic  force,  x de- 
notes the  motion  of  the  armature  in  the  direction  of  the 
magnetic  force,  Rt,  Lt,  Bc,  and  l denote,  the  resistance 
of  the  solenoid  coil,  the  inductance  of  the  coil,  the  mag- 
netic flux  density,  and  the  length  of  the  armature  which 
moves  orthogonal  to  the  magnetic  field,  respectively. 
m.  b,  and  k represent  the  effective  mass,  damping,  and 
the  stiffness  of  the  armature/poppet  system. 


The  fluidic  part  can  bc  modeled  using  the  unsteady 
Bernoulli  equation,  and  the  conservation  of  mass  across 
the  injector  assuming  incompressible  flow.  The  un- 
steady velocity,  v,  can  be  obtained  from  the  former  ap- 
plied between  the  inlet  to  the  valve  which  is  connected 
to  a pressurized  tank  (where  the  flow  velocity  is  « 0, 
and  p — p„),  and  the  outlet  to  the  combustor  where 
p = pc,  and  A p — p0  — pc,  as 

+ \l,u~  = (25) 

In  case  small  perturbations  in  pc  affect  the  velocity  out 
of  the  injector,  the  unsteady  velocity  out  of  the  valve,  v, 
is  linearized  as 


dv'  , 
Tfluid~dt  + V 


-Pc 

\/2Epp ' 


(26) 


where  t fluid  — Li  / \/2 A p / p is  the  fluidic  time  con- 
stant, and  Li  is  the  distance  between  the  tank  and  the 


valve’s  outlet.  TfiUid  is  negligible  for  conditions  where 
Lf  « lm,  A f>  is  large  (which  is  expected),  and  p is 
small  (0(iky/m3)  for  most  gaseous  fuels)6. 

The  mass  flow  rate,  defined  as  rhj  — pvA , is  perturbed, 
assuming  oscillations  in  v (caused  by  the  dynamics  in 
Eq.  (26))  and  ,1  (caused  by  the  motion  of  the  poppet  x), 
as 

m'f  — pAv'  — pvA' , (27) 

w'hcrc  wc  assume  A'  — k0x,  and  ka  > 0.  Equations 
(26)  and  (27)  describe  the  fluid  dynamics  due  to  per- 
turbations in  p{.  and/or  x.  In  most  practical  cases,  Aj7 
is  large  to  guarantee  choked  conditions  in  the  injector’s 
discharge,  thus,  the  first  term  in  the  RHS  of  Eq.  (27) 

can  be  neglected.  Using  v — from  Eq.  (25),  we 

simplify  Eq.  (27)  as 

m'f  “ k„pJ~-x.  (28) 


Equations  (21)-(24)  and  (28)  determine  the  input- 
output  relation  between  the  fuel-injector  input  E and 
the  output  m'f  which  is  expressed  in  the  Laplace  do- 
main as: 

ThW  = k no 

E{s)  (r€s  + -f-  bs  •+*/?)  + B~l2/Rt  $ 

where  kv  = Belk0p\J=^/Re. 


In  most  solenoid  systems,  the  armature  electric  time 
constant.  rt  = Lt./Rt,  is  negligible  compared  to  the 
acoustics  time  constant  [27].  In  these  valves,  the  stiff- 
ness of  the  spring,  k , is  large,  for  a fast  closing  of  the 
valve,  when  the  voltage  is  turned  off.  Also,  the  mass, 
to,  of  the  armature  is  very  small  in  many  of  the  typical 
injectors  to  minimize  inertia  forces  [27].  The  damping 
term,  b,  contains  the  overall  damping  including  stiction 
and  friction,  and  typically  is  large.  Thus,  the  mechan- 
ical system  can  bc  simplified  as  a first-order  system;  a 
damper-spring  system  [28].  The  mechanical  time  con- 
stant usually  limits  the  bandwidth  of  typical  injectors  to 
approximately  100  Hz.  (Note  that  other  effects,  such  as 
impact  dynamics  are  not  included  here,  since  we  expect 
them  to  be  of  higher  frequencies  than  the  combustor  dy- 
namics). Thus,  Eq.  (29)  is  simplified  as 


m'j  is)  kvrm 

E{s)  7~m-&  A 1 ’ 

where  rm  — (b  + B‘jl2 / Rc) / k 7 


(30) 


6 In  the  case  of  liquid  fuels,  the  time  constant  can  be  comparable 
to  the  acoustics  time  constant  due  to  large  p,  (J[  1 OOOkp/frk ) 

7Fot  more  advanced  proportional  injectors,  internal  feedback 
loops  exist  (using  for  example  a position  transducer  for  the  arma- 
ture) to  guarantee  accurate  metering,  and  increase  its  bandwidth,  e.g., 
a Moog  DDV  proportional  valve  has  a bandwidth  of  450  Hz  [29]. 
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Figure  2:  Block  diagram  of  a typical  two-position  injector. 


2.5.2  A Two-Position  (on-off)  Fuel  injector: 

Some  fuel  injectors  currently  used  for  combustion  con- 
trol [6,  11,  30]  operate  only  between  two  positions,  on 
and  off.  Unlike  the  proportional  injectors  discussed 
above,  the  physical  stops  play  a more  prominent  role 
in  the  dynamics.  However,  one  can  still  model  two- 
position  injectors  in  the  same  manner  as  above  by  in- 
cluding the  effect  of  the  physical  stops  as  a saturation 
block  together  with  Eq.  (30)  as  shown  in  Fig.  2.  A two 
position  injector  is  set  on , after  the  voltage  input  over- 
comes a certain  threshold,  thus  creating  a dead-zone  in 
the  control  input  (see  figure)  which  will  be  discussed 
further  when  control  is  implemented  in  Sec.  3. 

An  additional  point  to  note  is  the  distinction  between 
the  injector  dynamics  during  transition  from  closing 
and  opening.  Typically,  the  injector  is  over-driven  by 
a high  voltage  in  the  opening  mode  to  ensure  fast  open- 
ing. The  opening  time  constant  is  different  than  the 
closing  one.  This  effect  can  also  be  included  in  the  in- 
jector model,  as  seen  in  Fig.  2,  by  assuming  that  rm 
varies  between  two  values  rTOl  and  Tm.,  depending  on 
whether  the  injector  is  transitioning  from  off  to  on  or 
on  to  off.  Note  that  rmi  — rm  (as  defined  before)  and 
rm„  — b/k.  We  refer  the  reader  to  [31]  for  more  de- 
tails regarding  this  model  when  validated  against  tw'o 
different  Parker-Hannafin  Series-9  injectors  9-130-905 
and  9-633-900. s 


3 Control 

In  this  section,  we  investigate  model-based  control 
strategies  for  abating  combustion  instability  using  sec- 

8It  should  be  noted  that  when  inertia  forces  in  the  armature  are 
important,  a second-order  fuel-injector  model,  with  high  damping, 
is  more  appropriate.  The  experimental  measurements  in  [30]  for  a 
General  Valve  injector  show  similar  dynamics. 


ondary  fuel  injection.  We  assume  that  the  pressure  sig- 
nal is  the  measured  output,  using  a pressure  transducer. 
The  transducer  dynamics  is  neglected  since  it  typically 
has  a much  higher  bandw'idth  than  the  combustion  dy- 
namics. We  examine  the  control  with  an  injector  lo- 
cated at  (i)  the  burning  zone  or  (ii)  further  upstream. 
We  assume  in  particular  that  the  combustion  dynamics 
is  determined  by  several  coupled  acoustic  modes  [20], 
We  also  assume  that  instability  is  primarily  induced  by 
fluctuations  in  the  flame  area  coupled  with  the  acous- 
tics. We  refer  the  reader  to  [18, 19]  for  investigations  of 
instability  caused  by  equivalence  ratio  fluctuations  and 
its  control. 

3.1  Actuated  Combustor 

Denoting  the  contribution  of  the  fuel  injector  to  the 

equivalence  ratio  as  we  have  that  <p{.  — / a0, 

where  mc  and  ac,  are  the  mean  air  mass  flow  rate,  and 
the  stoichiometric  fuel  to  air  ratio,  respectively.  We  as- 
sume that  <p[.  is  uniform  radially,  and  that  perturbations 
are  carried  intact  by  the  mean  flow  to  the  burning  zone, 
after  a lime  delay  tc , where  r{,  = Lc/u.  Lc  is  the  dis- 
tance between  the  injector  discharge  and  the  burning 
zone,  and  u is  the  mean  velocity  of  the  reactants  in  the 
combustor. 

From  Eq.  (18),  the  impact  of  <p'.  on  the  combustion 
dynamics  can  be  taken  into  account  as 

iji  + 2C u>;  r/i  ~ ^kRu'  + chffft  - r<,)]  , (31) 

for  i = 1,2,..,  n.  where  i denotes  the  mode  number. 

3.2  Injection  at  the  Burning  Zone 

Injection  at  the  flame  has  shown  success  in  several  ex- 
perimental facilities  [6,  11,  30],  and  in  a practical  full- 
scale  170  MW  gas-turbine  combustor  [5].  We  carry 
out  active  control  design  assuming  that  the  injection  is 
at  the  flame,  using  the  model  in  Eq.  (31)  while  rc  ~ 0, 
together  with  the  injection  dynamics  described  by  Eq. 
(30).  The  input-output  model  relation  between  the  in- 
jector input  voltage,  E,  and  the  pressure,  p\  is  given 
by 


p'(s)  = W„(S)£(s),  (32) 

where  Ujft.s)  is  the  transfer  function  of  a finite- 
dimensional model  of  the  combustor.  kp,  Zp(s)  and 
Rp(s ) are  the  corresponding  gain,  numerator  and  de- 
nominator, respectively. 

An  optimal  control  for  the  finite-dimensional  system 
as  in  Eq.  (32)  is  LQG/LTR  [32].  This  method  has 
been  successful  in  suppressing  combustion  instability 
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[13].  Its  success  lies  in  its  ability  to  generate  sat- 
isfactory performance  over  a wide  range  of  frequen- 
cies. unlike  phase-shift  controllers  which  can  destabi- 
lize stable  dynamics  [33],  Experimental  validation  of 
the  LQG/LTR  for  combustion  control  has  been  demon- 
strated in  [31,  34,  35].  In  [34],  a similar  physically- 
based  model  was  used  with  a loudspeaker  as  an  actu- 
ator. In  [31,  35],  a system-ID  approach  based  oil  sub- 
space and  ARMAX  methods  [36]  was  used  to  suppress 
pressure  oscillations  in  a dump  and  a swirl-stabilized 
combustors,  respectively,  using  pulsed  injection.  In  this 
paper,  we  show  through  simulation  studies  that,  as  in 
[13],  LQG/LTR  can  also  be  used  successfully  based  on 
a physical  model,  using  fuel  injector  as  an  actuator.  For 
details  of  the  LQG/LTR  control  design,  we  refer  the 
reader  to  [13,  34]. 

3.2.1  Simulations  of  the  LQG/LTR  Controller: 

A fifth  order  combustor  dynamics  model  including  the 
first  two  modes,  the  flame  dynamics,  and  the  injector 
dynamics  is  considered.  The  combustor  parameters  and 
conditions  are  taken  as  in  [13],  these  cause  a three- 
quarter-mode  instability  which  resonates  at  500  Hz  ap- 
proximately, and  has  unsteady  pressure  amplitudes  of 
O(lOOPa). 

We  choose  first  a proportional  injector,  as  in  [5],  with 
a bandwidth  of  300  Ilz  which  is  in  the  range  of  avail- 
able high-speed  injectors  [29].  Figure  3 shows  the  time 
response  of  the  pressure  and  the  control  input,  </>'./© 
(<p  — 0.7).  Control  is  applied  at  t = 50 ms.  We  note 
that  although  the  bandwidth  is  lower  than  the  system 
dynamics,  the  control  with  proportional  injector  is  still 
capable  of  stabilizing  the  system,  since  the  dynamics  of 
the  injector  are  taken  into  consideration  in  the  design  of 
the  LQG/LTR.  As  seen  in  Eq.  (30).  the  injector  dynam- 
ics are  first  order,  and  the  gain  and  phase  introduced 
by  the  injector  depend  on  its  bandwidth.  For  smaller 
bandwidths  cases,  the  phase  increases  while  the  con- 
trol authority,  i.e.,  the  gain,  is  dramatically  decreased 
around  acoustic  frequencies.  The  controller,  however, 
has  enough  degrees  of  freedom  to  adjust  the  phase,  and 
increase  the  voltage  amplitude  into  the  injector  to  pro- 
duce the  required  control  authority,  around  the 

acoustics  frequencies  thus  guaranteeing  stability. 

While  a low  bandwidth  proportional  injector  is  capable 
of  maintaining  the  system  at  vanishingly  small  pressure 
perturbations,  a two-position  injector  may  not  be  as  ef- 
fective [11,  30,  35].  Using  a 300  Hz  bandwidth  injec- 
tor, the  LQG/LTR  is  still  capable  of  stabilizing  the  sys- 
tem, but  the  pressure  is  suppressed  to  a small  but  finite 
amplitude  limit  cycle  as  seen  in  Fig.  4.  The  reason 
is  that  the  injector  has  a threshold  input  voltage  value 
at  which  it  is  activated.  Thus,  following  the  suppres- 
sion of  the  instability,  the  injector  stops  pulsing  at  from 


Figure  3:  Response  of  the  controlled  combustor  with  a pro- 
portional injector  with  a bandwidth  of  « 300  Hz. 


Figure  4:  Response  of  the  controlled  combustor  with  an  on- 
off  injector  set  to  deliver  « 0.125  and 

with  a bandwidth  « 300  Hz. 


t «105-138ms,  as  seen  in  the  figure.  Disturbances  in 
the  combustor  force  the  pressure  to  grow',  until  the  mea- 
sured voltage  by  the  microphone  reaches  the  threshold 
at  which  the  injector  starts  to  fire  again.  In  the  case 
simulated,  this  occurs  at  t > 13 8ms.  This  sequence  is 
repeated  indefinitely. 

Increasing  the  control  fuel-flow'  rate,  and  thus  ©c,|mnj, 
the  combustor  is  stabilized  in  a smaller  settling  time.  As 
seen  in  Fig.  5,  when  4>'c\mux / <P  is  doubled,  the  settling 
time  diminishes  by  ~ 80%.  Moreover,  the  rms  of  the 
steady-state  pressure  is  smaller.  A similar  effect  has 
been  observed  in  [1 1], 

We  also  investigate  the  effect  of  bandwidth  for  a two- 
position  injector.  Limiting  it  to  50  Hz,  as  seen  in  Fig. 
6,  the  pressure  settles  to  a higher-amplitude  limit  cycle, 
and  the  injector  is  incapable  of  tracking  the  command 


Figure  5:  Response  of  the  controlled  combustor  with  on-off 
injector  set  to  deliver  isiszz.  « 0.27,  and  with  a 
bandwidth  of  300  Hz. 
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Figure  6:  Response  of  the  controlled  combustor  with  on- 
off  injector  with  lower  bandwidth  « 50  Hz  (the 
acoustics  unstable  frequency  is  500  Hz.  approxi- 
mately). 


from  the  controller;  the  injector  stays  open  all  the  time. 
This  show's  that  injector  bandwidth  is  a serious  prob- 
lem [11,  35j.  Different  solutions  have  been  proposed 
that  include:  (i)  Developing  faster  injectors  [29],  (ii) 
Use  of  multiple  injectors  w'hich  are  fired  alternatively 
to  increase  the  apparent  frequency  of  actuation  [11],  A 
different  approach  that  has  shown  promise  regardless 
of  high-bandwidth  injectors  is  through  fuel  pulsing  at 
low  frequencies  (much  lower  than  the  acoustics).  This 
is  demonstrated  experimentally  in  [4,  37],  and  analyti- 
cally in  [38]. 

3.3  Injection  Upstream  the  Burning  Zone;  Delay  in 
the  Control  Input 

While  injecting  fuel  directly  on  the  flame  [5,  6,  1 1,  30] 
can  avoid  actuation  delays,  it  introduces  hot  spots  at  the 
flame  surface  thus  increasing  emissions.  In  addition,  if 
mixing  is  weak  at  the  injection  port,  we  run  the  danger 
of  creating  a secondary  diffusion  flame  which  can  be 
completely  decoupled  from  the  main  premixed  flame, 
and  hence  become  ineffective  in  suppressing  the  insta- 
bility9. 

In  this  section,  we  study  the  effects  of  pulsed-fuel  in- 
jection upstream  the  burning  zone.  This  has  been  uti- 
lized in  [3]  where  secondary'  injection  was  done  at  the 
primary  fuel  source.  In  [19],  we  presented  a Posi-Cast 
control  capable  of  working  W'ith  an  injector  located  at 
an  arbitrary  distance  upstream  the  flame.  In  that  case,  a 
bulk  mode  was  unstable.  Here,  we  extend  the  analysis 
of  the  Posi-Cast  control,  and  show  that  it  is  capable  of 
stabilizing  longitudinal  modes  as  well. 

3.3.1  Posi-Cast  Control:  A powerful  approach 

for  controlling  systems  with  known  time-delay  W'as 
originated  by  Smith  [40],  know'n  also  as  Posi-Cast  for 
“positive  forecasting"  of  future  states.  The  idea  is  to 
compensate  for  the  delayed  output  using  input  values 
stored  over  a time  w'indow'  equal  to  the  delay  time,  i.e. 
[f  — Tc,t ],  and  estimate  the  future  output  using  a model 

9Tliis  has  been  noticed  in  experiments  at  MIT  and  at  UTRC  [39], 


of  the  combustor.  Only  stable  systems  were  consid- 
ered. An  extension  to  include  unstable  systems  was 
proposed  in  [41]  using  finite-time  integrals  of  the  de- 
layed input  values  thereby  avoiding  unstable  pole-zero 
cancellations  which  may  occur.  A frequency-domain 
pole-placement  technique  for  unstable  systems  was  first 
proposed  in  [42]  and  a similar  technique  will  be  pre- 
sented here. 

The  model  in  Eq.  (32),  in  the  presence  of  a time  delay, 
r,.,  can  be  re-written  as 

p'(t)  = Wp{s){E(t  - n,)J,  n>(«)  - ^^33) 


Due  to  the  nature  of  the  combustion  system,  not  all 
states  are  accessible,  only  the  system  input,  i.e.  the  volt- 
age to  the  injector  17(f),  and  the  output,  p'  in  our  case 
are  measured.  A standard  pole-placement  controller  is 
required  (for  more  information,  see  [42, 43]).  The  pres- 
ence of  the  time-delay,  rc.  in  the  control  input,  moti- 
vates the  use  of  an  additional  signal  in  the  control  in- 
put, E(t),  denoted  as  Ei  (t)  which  anticipates  the  future 
output  using  a model  of  the  system  [19],  The  resulting 
controller  structure  is  described  as 


m 

£i(f) 


(34) 


TJ, 


Rp  («) 


Rp(s) 


where  A (s)  is  a chosen  stable  polynomial  of  degree  «•  — 
1,  d(s) , n-i  (s)  and  n 2 (s),  are  polynomials  of  degree  n — 
1 at  most,  and  c(s)  is  of  degree  n — 2 at  most.  For 
stability,  these  must  satisfy  the  relations 


c($)Rp(s)  + kpd[s)Zp{s)  — A(s)n-2(s),  (35) 

n-i(s)  = Rp(-s)  - Rm.(s $36) 

where  Rm(s ) is  the  desired  characteristic  equation, 
which  is  a stable  monic  polynomial  of  the  same  order 
of  Rp(s). 

Using  the  controller  structure  in  Eq.  (34)  with  the  con- 
ditions in  Eqs.  (35)  and  (36),  the  closed-loop  transfer 
function  can  be  computed  as 


Wcl{s) 


kpe  src 
Rm(s)  ' 


(37) 


The  control  input  law.  in  Eq.  (34),  introduces  additional 
dynamics  including  non-minimum  phase  zeros  having 
the  same  eigen-values  of  Rp(s).  Obviously,  these  lead 
to  unstable  pole-zero  cancellations  since  the  combus- 
tor model  is  open-loop  unstable  (i.e.,  Rp(s ) has  unsta- 
ble eigen  values).  Unstable  pole-zero  cancellations  are 


Figure  7:  Response  of  the  controlled  combustor  with  a time- 
delay  of  100ms  in  the  input  signal,  proportional 
injector.  Note:  only  the  envelope  of  the  response 
is  shown  for  clarity,  since  the  scale  of  the  plot  does 
not  permit  seeing  individual  cycles. 


known  to  cause  problems  concerning  observability  and 
controllability  of  the  plant  (see  [27]  for  more  details). 
As  a result,  a modification  in  the  synthesis  of  Ei(t)  in 
Eq.  (34)  was  suggested  by  [41],  To  avoid  unstable  pole- 
zero  cancellations,  £)  (t)  must  be  generated  as  a finite 
integral  of  the  form 

Ex (t)  = E(/0  e-*'J?(t  + a)*r),(38> 

where  A,  ’s  are  the  eigen  values  of  the  combustor  sys- 
tem, i.e.  Rp(. s)  = n]L1(.s  — A;).  Taking  the  Laplace 
transform  of  Eq.  (38  ),  one  can  show  that 

nx(s)  _ -yp  <*»  _ A flj 

Rp(tt)  “ » — A,:  ’ Rp{a)  “ s — A,:  ’ 

w here  /7,  = cv,:eXirc.  Another  condition  for  the  success- 
ful use  of  the  finite  integral  in  Eq.  (38)  is  that  Rp(s)  has 
no  repeated  roots  [42], 

The  controller  described  in  Eqs.  (34)  and  (38)  is  suf- 
ficient to  stabilize  the  combustor  provided  that  an  ac- 
curate description  of  the  plant  and  the  time  delay  are 
available.  This  controller  has  been  show'll  to  provide 
robustness  to  uncertainties  in  the  plant  including  the 
time  delay  [41].  Adaptive  versions  of  the  same  con- 
troller have  been  investigated  [44,  45],  and  have  shown 
to  extend  the  robustness  of  the  controller  to  parameter 
uncertainties. 

3.3.2  Simulations  of  the  Posi-Cast  Controller: 

The  controller  in  Eqs.  (34)  and  (38)  is  implemented  for 
injection  at  a distance  of  ~ 3cm  upstream  the  burning 
zone.  tc  is  estimated  to  be  1 00ms,  which  is  about  50 
times  the  time  constant  of  the  unstable  frequency. 

The  closed-loop  simulation  is  illustrated  in  Fig.  7.  Al- 
though control  is  switched  on  at  t — 50 ms,  the  pressure 
keeps  increasing  for  an  additional  /,  = tc  — 100ms 
(from  t — 50—  150ms),  then  stalls  for  another  1 00  ms 
(from  t — 150  — 250ms)  before  decaying.  The  reason 


Figure  8:  Response  of  the  controlled  combustor  with  a time- 
delay  of  100  ms  in  the  input  signal,  on-off  injector. 
Note:  only  the  envelope  of  the  response  is  shown 
for  clarity',  since  the  scale  of  the  plot  does  not  per- 
mit seeing  individual  cycles. 


for  the  former  delay  is  physical  and  is  due  to  the  time 
taken  for  the  pulsed-fuel  to  reach  the  burning  zone.  The 
latter  is  due  to  a computational  delay  in  the  controller. 
Specifically,  the  finite  integral  in  Eq.  (38)  outputs  in- 
correct values  for  a period  of  tv.  This  is  because  the 
computation  of  the  finite  integral  relies  on  a stored  win- 
dow of  the  past  values  of  the  control  input  of  the  size  of 
rr.  When  control  is  switched  on,  the  window  consists 
of  control  inputs  proportional  to  p'  which  has  not  yet 
“felt”  the  effect  of  control  due  to  the  physical  delay  rc 
(the  values  of  p'  are  still  those  of  the  open-loop  com- 
bustor). It  requires  therefore  t — 2 tc  to  start  forming 
a window  of  integration  with  control  input  correspond- 
ing to  closed-loop  values.  This  confirms  observations 
in  [41]. 

In  Fig.  8,  a two-position  injector  is  used.  The  control 
design  is  based  on  the  linear  model,  and  its  parameters 
are  fine-tuned  to  handle  the  nonlinearities.  As  discussed 
earlier,  the  control  is  switched  on  at  50ms,  and  stabi- 
lizes the  system.  The  injector  stays  on  as  long  as  the 
voltage  signal  into  the  injector  is  greater  than  a thresh- 
old. as  discussed  before  in  Sec.  2.5.2. 

It  should  be  noted  that  when  combustion  instability  is 
caused  by  fluctuations,  as  discussed  in  Sec.  2,  the 
characteristic  equation  will  look  different  than  in  Eqs. 
(32)  and  (33).  Rp(s)  will  have  terms  which  are  delayed, 
due  to  the  convective  delay',  r5,  carried  by  Hence, 
Rp{s)  becomes  infinite  dimensional.  To  circumvent 
this,  a Pade  approximation  [46]  is  used  to  get  a finite  di- 
mensional description  of  Rp{s).  and  thus  the  LQG/LTR 
and  Posi-Cast  controllers  as  described  in  Secs.  3.2  and 
3.3.  respectively,  can  similarly  be  used  for  this  case. 


4 Summary 

In  this  paper,  a complete  model  of  the  combustion  dy- 
namics leading  to  instability  is  developed.  A model 
encompassing  the  acoustics,  the  heat  release,  the  cou- 
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pling,  mixing  and  injector  dynamics  is  presented.  The 
heat  release  is  derived  using  flame  kinematics  for 
flames  with  high  Damkohler  numbers  and  moderate  tur- 
bulence, and  the  effect  of  forcing  in  the  velocity  and  the 
equivalence  ratio  is  illustrated.  Different  stability  cri- 
teria are  discussed  for  the  different  dynamics  resulting 
from  the  coupling  of  these  with  the  acoustics.  While 
velocity  perturbations  cause  instability  of  a phase-lag 
nature,  equivalence-ratio  oscillations  introduce  a time- 
delay  instability.  Injection  both  at  and  upstream  of  the 
burning  zone  is  studied.  An  LQG/LTR  control  is  used 
for  the  former  while  a Posi-Cast  control  is  developed 
and  used  for  the  latter.  Both  proportional  and  two- 
position  (on-off)  pulsed  injection  is  examined.  Pro- 
portional injectors  show  superior  performance  to  two- 
position  ones.  Performance  is  shown  to  be  plagued  by 
limited  injector  authority  and  bandwidth.  Some  of  the 
short-term  fixes  proposed  for  this  problem  include  the 
use  of  multiple  injectors  [ 1 1 ] and  low  frequency  pulsing 
[4,  37,  38],  A long-term  solution  is  clearly  the  design 
of  a high-speed,  high-authority  injector. 
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PAPER  -23,  J.  P.  Hathout 

Question  (A.  Banaszuk.  USAt 

How  effective  is  the  Smith  predictor-type  controller  if  there  is  a large  unmeasurable  dis- 
turbance, such  as  turbulent  velocity  components  driving  the  system?  Don’t  you  use  the 
assumption  that  you  know  all  the  inputs  during  the  delay  block? 

Reply 

Robustness  of  Posi-Cast  controllers  has  been  analyzed  in  several  publications.  See,  for 
example,  references  [41]  and  [45]  in  the  paper. 

Question  ( J.  Hermann.  Germany ) 

As  you  have  shown,  in  a system  with  long  time  delay,  there  must  be  a long  delay  from 
activating  the  active  control  to  achieving  suppression  of  the  oscillations.  Applied  to  a 
realistic  system,  this  may  be  too  long.  Can  this  problem  be  reduced? 

Renlv 

The  large  delay  vc  is  due  to  the  fact  that  the  secondary  injection  is  at  a certain  distance 

upstream  of  the  burning  zone.  This  is  a physical  constraint  that  affects  closed-loop  con- 
trol performance  in  general,  regardless  of  the  control  algorithm.  The  implication  of  this 
constraint  is  that  the  effect  of  the  control  on  the  pressure  will  not  take  place  for  a mini- 
mum period  of  rc . After  this  period,  the  more  sophisticated  the  controller,  the  smaller  the 
settling  time.  The  Posi-Cast  control  indeed  results  in  the  latter,  as  described  in  the  pres- 
entation. 

The  time-delay  effect  may  be  reduced  by  introducing  other  actuating  inputs  closer  to  the 
burning  zone.  For  example,  it  may  be  possible  to  introduce  acoustic  drivers  or  auxiliary 
fuel-injectors  closer  to  the  flame,  which  will  not  have  the  same  time-delay  constraint. 

The  disadvantage  of  injection  close  to  the  flame  may  lie  in  the  increase  of  emissions;  an 
acoustic  driver  obviously  avoids  this  problem,  but  may  have  limited  authority.  In  either 
case,  by  designing  a multi-input  control  strategy,  where  inputs  closer  to  the  flame  act 
only  over  the  time  interval  [0,TC  ],  it  may  be  possible  to  cope  with  the  large  delay. 


